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Abstract - We introduce a quantum version of the Game of Life and we use it to study the emer- 
gence of complexity in a quantum world. We show that the quantum evolution displays signatures 
of complex behaviour similar to the classical one, however a regime exists, where the quantum 
Game of Life creates more complexity, in terms of diversity, with respect to the corresponding 
classical reversible one. 



The Game of Life (GoL) has been proposed by Conway 
in 1970 as a wonderful mathematical game which can de- 
scribe the appearance of complexity and the evolution of 
"life" under some simple rules pQ. Since its introduction 
it has attracted a lot of attention, as despite its simplicity, 
it can reveal complex patterns with unpredictable evolu- 
tion: From the very beginning a lot of structures have 
been identified, from simple blinking patterns to complex 
evolving figures which have been named "blinkers" , "glid- 
ers" up to "spaceships" due to their appearance and/or 
dynamics j2]. The classical GoL has been the subject of 
many studies: It has been shown that cellular automata 
defined by the GoL have the power of a Universal Turing 
machine, that is, anything that can be computed algo- 
rithmically can be computed within Conway's GoL [3111]. 
Statistical analysis and analytical descriptions of the GoL 
have been performed; many generalisations or modifica- 
tions of the initial game have been introduced as, for exam- 
ple, a simplified one dimensional version of the GoL and a 
semi-quantum version [5H7J . Finally, to allow a statistical 
mechanics description of the GoL, stochastic components 
have been added [8]. 

In this letter, we bridge the field of complex systems 
with quantum mechanics introducing a purely quantum 
GoL and we investigate its dynamical properties. We show 
that it displays interesting features in common with its 
classical counterpart, in particular regarding the variety of 
supported dynamics and different behaviour. The system 
converges to a quasi-stationary configuration in terms of 
macroscopic variables, and these stable configurations de- 
pend on the initial state, e.g. the initial density of "alive" 
sites for random initial configurations. We show that sim- 
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Fig. 1: Example of the evolution of the GoL described by 
Hamiltonian JT]) for a simple initial configuration. Empty 
(blue) squares are "dead" sites, coloured (red) ones are "alive" . 



pie, local rules support complex behaviour and that the 
diversity of the structures formed in the steady state re- 
sembles that of the classical GoL, however a regime exists 
where quantum dynamics allows more diversity to be cre- 
ated than possibly reached by the classical one. 

The universe of the original GoL is an infinite two- 
dimensional orthogonal grid of square cells with coordi- 
nation number eight, each of them in one of two possible 
states, alive or dead pQ. At each step in time, the pattern 
present on the grid evolves instantaneously following sim- 
ple rules: any dead cell with exactly three live neighbours 
comes to life; any live cell with less than two or more than 
three live neighbours dies as if by loneliness or overcrowd- 
ing. As already pointed out in [7J, the rules of the GoL are 
irreversible, thus their generalisation to the quantum case 
implies rephrasing them to make them compatible with a 
quantum reversible evolution. The system under study is 
a collection of two-level quantum systems, with two pos- 
sible orthogonal states, namely the state "dead" (|0)) and 
"alive" (|1)). Clearly, differently from the classical 
site can be also in a superposition of the two possible clas- 
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Fig. 2: Colour on-line. From left to right: Countour plot of the time evolution of the populations (ni(t)) (column 1), vis- 
ibility Vi(t) (column 2), discretized populations T>i(t) (column 3) and clustering C(£,t) (column 4) for three different initial 
configurations: four alive sites separated by two dead ones (A), twenty- four alive sites grouped together (B) and a random 
initial configuration (C). Time is reported on the x-axis (in arbitrary units), and position (cluster size) i = 1, ...,L on the 
y-axis in columns one to three (four). Arrows in panel 3A highlight the three subsequent generations of a "blinker" reported 
schematically in Fig. [3] The colour code goes from zero to M = 1 (M = 4 for the clustering and to M = .1 for the visibility), 
from blue through green to red. 



sical states. The dynamics is denned as follows in terms of 
the GoL language: a site with two or three neighbouring 
alive sites is active, where active means that it will come 
to life and eventually die on a typical timescale T (setting 
the problem timescale, or time between subsequent gener- 
ations). That is, if maintained active by the surrounding 
conditions, the site will complete a full rotation, if not, it 
is "frozen" in its state. Stretching the analogy with Con- 
way's GoL to the limit, we are describing the evolution of 
a Virus culture: each individual undergoes its life cycle if 
the environment allows it, otherwise it hibernates in its 
current state and waits for conditions to change such that 
the site may become active again. This slight modification 
allows us to recover the reversibility of the dynamics and 
to introduce a quantum model that, as we shall see, re- 
produces most of the interesting complex behaviour of the 



classical GoL from the point of view of a classical observer. 
However, its evolution is purely quantum and thus we are 
introducing a tool that will allow to study the emergence 
of complexity from the quantum world. 

Model. — The Hamiltonian describing the aforemen- 
tioned model is given by 

# = X> + ^)-(A/?+A/?) (1) 

f=3 

where L is the number of sites; b and are the usual an- 
nihilation and creation operators (h = 1); the operators 
Mf = J2p n a npn 7 fis and = J2p> n a npn-yfis (n = tfb, 
n = 1 — n, the indices a, /3, 7, 6 label the four neighbouring 
sites) count the population present in the four neighbour- 
ing sites (the sum runs on every possible permutation P 
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Fig. 3: Schematic representation of a one-dimensional time- 
evolution of the discretized population T>i{t) of a "blinker" 
(case A of Fig. [2}. From left to right the states of subsequent 
generations are sketched. Empty (blue) squares are "dead" 
sites, coloured (red) ones are "alive" 



and P 1 of the positions of the n and n operators) and N 2 
(A/" 3 ) gives the null operator if the population is differ- 
ent from two (three), the identity otherwise. For classical 
states, as for example an initial random configuration of 
dead and alive states, the Hamiltonian ([1]) is, at time zero, 
H Active = bi + b\ on the sites with two or three alive neigh- 
bours and Hnibemate = otherwise. If the Hamiltonian 
would remain constant, every active site would oscillate 
forever while the hibernated ones would stand still. On the 
contrary as soon as the evolution starts, the state evolves 
into a superposition of possible classical configurations, 
resulting in a complex dynamics as shown below and the 
interaction between sites starts to play a role. Thus, the 
Hamiltonian introduced in Eq. ([lj induces a quantum dy- 
namics that resembles the rules of the GoL: a site with 
less than two or more than three alive neighbouring sites 
"freezes" while, on the contrary, it "lives". The difference 
with the classical game - connected to the reversibility 
of quantum dynamics - is that "living" means oscillat- 
ing with a typical timcscale between two possible classical 
states (see e.g. Fig. [TJ 

Dynamics. — To study the quantum GoL dynamics 
we employ the time dependent Density Matrix Renor- 
malization group (DMRG). Originally developed to in- 
vestigate condensed matter systems, the DMRG and its 
time dependent extension have been proven to be a very 
powerful method to numerically investigate many-body 
quantum systems [9HT2]. As it is possible to use it effi- 
ciently only in one-dimensional systems, we concentrate to 
the one-dimensional version of the Hamiltonian ([T]): the 
operators A/" 2 and A/" 3 count the populated sites on the 
nearest-neighbour and next-nearest-neighbour sites and 
thus a = i — 2,(3 = i — 1,7 = i + 1,6 = i + 2. Note 
that it has been shown that the main statistical proper- 
ties of the classical GoL are the same in both two- and 
one-dimensional versions [6]. 

To describe the system dynamics we introduce different 
quantities that characterise in some detail the system evo- 
lution. We first concentrate on the population dynamics, 
measuring the expectation values of the number operator 
at every site (rn(t)). This clearly gives a picture of the 
"alive" and "dead" sites as a function of time, as it gives 
the probability of finding a site in a given state when mea- 




Fig. 4: Left: Average population p(t) (upper) and diversity 
A(f) (lower) as a function of time for different initial popula- 
tion density po- Right: Equilibrium average population p (up- 
per) and diversity A (lower) for the quantum (blue squares) 
and classical (red circles) GoL as a function of the initial pop- 
ulation density po- Simulations are performed with a t-DMRG 
at third order, Trotter step St = 10~ 2 , truncation dimension 
m = 30, size L = 32, averaged over up to thirty different initial 
configurations. 



sured. That is, if we observe the system at some final time 
Tf we will find dead or alive sites according to these prob- 
abilities. In Fig. [2] we show three typical evolutions (left- 
most pictures): configuration A corresponds to a "blinker" 
where two couples of nearest-neighbour sites oscillate reg- 
ularly between dead and alive states (a schematic rep- 
resentation of the resulting dynamics of the discretised 
population T>i (t) is reproduced also in Fig. [3]) ; configura- 
tion B is a typical overcrowded scenario where twenty-four 
"alive" sites disappear leaving behind only some residual 
activity; finally a typical initial random configuration (C) 
is shown. Notice that in all configurations it is possible 
to identify the behaviour of the wave function tails that 
propagate and generate interference effects. These effects 
can be highlighted by computing the visibility of the dy- 
namics, the maximum variation of the populations within 
subsequent generations, defined as: 

T T 

Vi{t) = I max 71,(0 - min ni(t') |; * € [*- + ^ ^ 

that is, the visibility at time t reports the maximum vari- 
ation of the population in the time interval of length T 
centered around t. The visibility clearly follows the pre- 
ceding dynamics (see Fig. [2j second column) and identifies 
the presence of "activity" in every site. 

To stress the connections and comparisons with the orig- 
inal GoL we introduce a classical figure of merit (shown in 
the third column of Fig. [2]): we report a discretized ver- 
sion of the populations as a function of time (T*i(t) = 1 
for rii(t) > 0.5 and T>i(t) = otherwise). Notice that 
T>i(t) gives the most probable configuration of the system 



p-3 



D. Bleh T. Calarco S. Montangero 



after a measurement on every site in the basis {|0), |1)}. 
Thus, we recover a "classical" view of the quantum GoL 
with the usual definition of site status. For example, con- 
figuration A is a "blinker" that changes status at every 
generation (see FigJ3] and More complex configura- 
tions appear in the other two cases. The introduction of 
the discretized populations T>i can also be viewed as a 
new definition of "alive" and "dead" sites from which we 
could have started from the very beginning to introduce a 
stochastic component as done in [8] . This quantity allows 
analysis to be performed as usually done on the classical 
GoL and to stress the similarities between the quantum 
and the classical GoL. Following the literature to quan- 
tify such complexity, we compute the clustering function 
C(£,t) that gives the number of clusters of neighbouring 
"alive" sites of size I as a function of time [B]. For ex- 
ample, the function C(£,t) for a uniform distribution of 
"alive" sites would be simply C(L) = 1 and zero otherwise 
while a random pattern would result in a random cluster 
function. This function characterises the complexity of 
the evolving patterns, e.g. it is oscillating between zero- 
and two-size clusters for the initial condition A, while it is 
much more complex for the random configuration C (see 
Fig. HI rightmost column). 

Statistics. — To characterise the statistical proper- 
ties of the quantum GoL we study the time evolution of 
different initial random configurations as a function of the 
initial density of alive sites. We concentrate on two macro- 
scopic quantities: the density of the sites that if measured 
would with higher probability result in "alive" states 

p(t)=X>(t)/L; (3) 

i 

and the diversity 

A(f) = £C&t), (4) 

i 

the number of different cluster sizes that are present in the 
systems, that quantifies the complexity of the generated 
dynamics [5J|B]. Typical results, averaged over different 
initial configurations, are shown in Fig. @] (left). As it 
can be clearly seen the system equilibrates and the den- 
sity of states as well as the diversity reach a steady value. 
This resembles the typical behaviour of the classical GoL 
where any typical initial random configuration eventually 
equilibrates to a stable configuration. Moreover, we com- 
pare the quantum GoL with a classical reversible version 
of GoL corresponding to that introduced here: at every 
step a cell changes its status if and only if within the first 
four neighbouring cells only three or two are alive. Notice 
that, the evolution being unitary and thus reversible, the 
equilibrium state locally changes with time, however the 
macroscopic quantities reach their equilibrium values that 
depend non trivially only on the initial population den- 
sity. In fact, for the classical game, we were able to check 
that the final population density is independent of the 



system size while the final diversity scales as L 1 / 2 (up to 
2 10 sites, data not shown). Moreover, the time needed to 
reach equilibrium is almost independent of the system size 
and initial population density. These results on the scal- 
ing of classical system properties support the conjecture 
that our findings for the quantum case will hold in gen- 
eral, while performing the analysis for bigger system sizes 
is highly demanding. A detailed analysis of the size scal- 
ing of the system properties will be presented elsewhere. 
In Fig. [4] we report the final (equilibrium) population den- 
sity (right upper) and diversity (right lower) as a function 
of the initial population density for both the classical and 
the quantum GoL for systems of L = 32 cells. The equi- 
librium population density p is a non linear function of 
the initial one po in both cases: the classical one has an 
initial linear dependence up to half-filling where a plateau 
is present up to the final convergence to unit filling for 
po = L Indeed, the all-populated configuration is a stable 
system configuration. The quantum GoL follows a simi- 
lar behaviour, with a more complex pattern. Notice that 
here a first signature of quantum behaviour is present: the 
steady population density reached by the quantum GoL is 
always smaller than its classical counterpart. This is prob- 
ably due to the fact that the evolution is not completely 
captured by this classical quantity: the sites with popula- 
tion below half filling, i.e. the tails of the wave functions, 
are described as unpopulated by T>i. However, this missing 
population plays a role in the evolution: within the over- 
all superposition of basis states, a part of the probability 
density (corresponding to the states where the sites are 
populated) undergoes a different evolution than the clas- 
sical one. In general, the quantum system is effectively 
more populated than the classical p indicates. This differ- 
ence in the quantum and classical dynamics is even more 
evident in the dependence of the equilibrium diversity on 
the initial population density pq. In the classical case the 
maximum diversity is slightly above three: on average, in 
the steady state, there are no more than about three dif- 
ferent cluster sizes present in the system independently of 
the initial configuration. On the contrary -in the quan- 
tum case- the maximal diversity is about four, increasing 
the information content (the complexity) generated by the 
evolution by about 10 — 20%. 

These findings are a signature of the difference between 
quantum and classical GoL. In particular we have shown 
that the quantum GoL has a higher capacity of generating 
diversity than the corresponding classical one. This prop- 
erty arises from the possibility of having quantum super- 
positions of states of single sites. Whether purely quantum 
correlations (entanglement) play a crucial role is under in- 
vestigation. Similarly, as there is some arbitrariness in 
our definition of the quantum GoL, the investigation of 
possible variations is left for future work. 

The investigation presented here fits perfectly as a sub- 
ject of study for quantum simulators, like for example cold 
atoms in optical lattices. Indeed, the five-body Hamilto- 
nian (p} can be written in pseudo spin-one-half operators 
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(Pauli matrices) and thus it can be simulated along the 
lines presented in |13j . In particular, these simulations 
would give access to investigations in two and three di- 
mensions that are not feasible by means of t-DMRG [TO] . 

In conclusion we note that this is one of the few available 
simulations of a many-body quantum game scalable in the 
number of sites [TUTTO] . With a straightforward generali- 
sation (adding more than one possible strategy defined in 
Eq. ([TJ) one could study also different many-player quan- 
tum games. This approach will allow different issues to 
be studied related to many-player quantum games such as 
the appearance of new equilibria and their thermodynam- 
ical properties. Moreover, the approach introduced here 
shows that one might investigate many different aspects of 
many-body quantum systems with the tools developed in 
the field of complexity and dynamical systems: In partic- 
ular, the relations with Hamiltonian quantum cellular au- 
tomata in one dimension and quantum games [TU[T7] . Fi- 
nally, the search for the possible existence of self-organised 
criticality in these systems along the lines of similar inves- 
tigations in the classical GoL [TH] , if successful, would be 
the first manifestation of such effect in a quantum sys- 
tem and might have intriguing implications in quantum 
gravity [19j[2Q]. 

After completing this work we became aware of another 
work on the same subject [2Tj . 
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